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Discovery of photon index saturation in the black hole binary 

GRS 1915+105 

Lev Titarchuk 1 and Elena Seifina 2 
ABSTRACT 

We present a study of the correlations between spectral, timing properties and 
mass accretion rate observed in X-rays from the Galactic Black Hole (BH) binary 
GRS 1915+105 during the transition between hard and soft states. We analyze all 
transition episodes from this source observed with Rossi X-ray Timing Explorer 
(RXTE), coordinated with Ryle Radio Telescope (RT) observations. We show 
that broad-band energy spectra of GRS 1915+105 during all these spectral states 
can be adequately presented by two Bulk Motion Comptonization (BMC) compo- 
nents: a hard component (BMC1, photon index Ti = 1.7 — 3.0) with turnover at 
high energies and soft thermal component (BMC2, T 2 = 2.7—4.2) with character- 
istic color temperature < 1 keV, and the redskewed iron line (LAOR) component. 
We also present observable correlations between the index and the normalization 
of the disk "seed" component. The use of "seed" disk normalization, which is 
presumably proportional to mass accretion rate in the disk, is crucial to establish 
the index saturation effect during the transition to the soft state. We discovered 
the photon index saturation of the soft and hard spectral components at values of 
t 4.2 and 3 respectively. We present a physical model which explains the index- 
seed photon normalization correlations. We argue that the index saturation effect 
of the hard component (BMC1) is due to the soft photon Comptonization in the 
converging inflow close to BH and that of soft component is due to matter accu- 
mulation in the transition layer when mass accretion rate increases. Furthermore 
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we demonstrate a strong correlation between equivalent width of the iron line 
and radio flux in GRS 1915+105. In addition to our spectral model components 
we also find a strong feature of "blackbody-like" bump which color temperature 
is about 4.5 keV in eight observations of the intermediate and soft states. We 
discuss a possible origin of this "blackbody-like" emission. 

Subject headings: accretion, accretion disks — black hole physics — stars: individual 
(GRS 1915+105) :radiation mechanisms: non-thermal — physical data and pro- 
cesses 



1. Introduction 

The study of the characteristic changes in spectral and variability properties of X- 
ray binaries is proved to be a valuable source of information on the physics governing the 
accretion processes and on the fundamental parameters of black holes (BHs). 

The simultaneous study of the spectral and timing evolution of a BH source during 
a state transition has been a subject of many investigations [see references in a review by 



Remillard & McClintock (2006)]. Fender & Belloni (2004), hereafter FB04, introduced a 



classification of the spectral states in GRS 1915 +105 and studied the spectral state evolution. 
Using X-ray colors (hardness ratio) they introduced three spectral states. State A: in which 
the strong blackbody-like (BB) component of color temperature > 1 keV dominates in the 
overall spectrum and little time variability is detected. State B: similar to state A but 
substantial red-noise variability on scales > 1 s occurs in this state. State C: the spectra are 
harder than those in states A and B. Photon indices of the power-law components vary from 
1.8 to 2.5. White-red noise (WRN) variability on scales > 1 s takes place in this state. 

Furthermore FB04 discussed the connection between states A, B, C observed in GRS 
1915 with the three "canonical" states in black hole candidates (BHCs) also identified by 
their timing and spectral properties. 

At a low luminosity state the energy spectrum is dominated by a hard Comptonization 
component combined (convolved) with a weak thermal component. The spectrum of this low 
(luminosity) hard state (LHS) is presumably a result of Comptonization (upscattering) of 
soft photons, originated in a relatively weak accretion disk, off electrons of the hot ambient 



plasma [see e.g. Sunyaev & Titarchuk (1980)]. Variability in LHS is high (fractional 



root-mean-square variability is up to 40%) and presented by a flat-top broken power law 
(white-red noise) shape, accompanied by quasi-periodic oscillations (QPOs) in the range of 
0.01-30 Hz, observed as narrow peaks in the power density spectrum (PDS). In high soft 
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state (HSS) a photon spectrum is characterized by a prominent thermal component which is 
probably a signature of a strong emission coming from a geometrically thin accretion disk. 
A weak power-law component is also present at the level of not more than 20% of the total 
source flux. In the HSS the flat-top variability ceases, QPOs disappear and PDS acquires a 
pure power-law shape. The total variability in HSS is usually about 5% fractional rms. The 
intermediate state (IS) is a transitional state between LHS and HSS. Note in addition to 
LHS, IS and HSS sometimes very soft state (VSS) is observed in which the BB component is 
dominant and the power-law component is either very weak or absent at all. The bolometric 
luminosity in VSS is a factor of 2-3 lower than that in HSS. 

FB04 concluded that probably all three states A, B, C of GRS 1915+105 are instances 
of something similar to the HSS/IS observed in other BHC systems, associated to the high 
accretion rate value for this source, although during the hardest intervals LHS might be 
sometimes reached. We come to the similar conclusions analyzing spectral and timing data 
from GRS 1915+105 obtained by RXTE (see below). 

Close correlations of nearly periodic variability [quasi-periodic oscillations (QPO)] ob- 
served during low-hard and intermediate states with the photon index of the Comptonization 
spectral component have been reported in multiple state transitions observed from accreting 



BHs [see |Vignarca et al. | ( [2003] ) , |Shaposhnikov fc Titarchuk] p006| ) , (2007), (2009), hereafter 
V03, ST06, ST07 ST09 respectively]. The ubiquitous nature of these correlations suggests 
that the underlying physical processes which lead to the observed variability properties are 
closely tied to the corona; furthermore, they vary in a well defined manner as the source 
makes a transition between spectral states. The fact that the same correlations are seen in 
many sources, which vary widely in both luminosity (presumably with mass accretion rate) 
and state, suggests that the physical conditions controlling the index and the low frequency 
QPOs are characteristics of these sources. Moreover, they may be an universal property 



of all accreting compact systems, including neutron sources too [see Titarchuk, & Fiorito 



(2004), hereafter TF04, and Titarchuk k Shaposhnikov (2005) 



When a BH is in LHS, radio emission is also detected and a jet is either seen or inferred 
( Fender et al. ||2001 ). Several models are successful in reproducing the energy spectrum from 



the radio domain to the hard X-rays [see e.g. Markoff et al. (2001 ), Vadawale et al. (2001 ), 



(2003 


) and 


Giannios 


(2005) 



The multiplicity of 

models that can fit well the time average spectrum of galactic BHs indicates that this alone 
is not enough to distinguish the most realistic one among them. X-ray timing features can 
be the key features to finally finding the common physical connection between the corona, 
the accretion disk and the jet radio emission in BHs. 



There is a big debate in the literature on the origin of quasiperiodic oscillation (QPO) 
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frequencies [see e.g. Remillard & McClintock (2006)] and its connection with the radio 
emission. Migliari et al. (2005) reported on correlation between radio luminosity and X-ray 



timing features in X-ray binaries containing a number of low magnetic field neutron stars and 
one black hole GX 339-4. They showed that in the low-hard state (LHS) radio luminosity 
is correlated with the low frequency QPO (LFQPO). Note ST09 demonstrate that in LHS 
of Galactic BHs LFQPO changes by order of magnitude, from 0.2 to 2 Hz whereas the 
photon index has almost the same value of 1.5. Below we show that in GRS 1915+105 the 
photon index monotonically increases with LFQPO and disk mass accretion rate, although 
the radio luminosity does not correlate with LFQPO and X-ray luminosity in the whole 
range of spectral states, from low-hard to high-soft through intermediate states. Recently 



Kylafis et al. (2008) suggested a model which explains how the QPO phenomenon is related 



to an appearance of radio flares (jets). Below (see §3) we present details of our observational 
study of the QPO connection with the X-ray and radio flaring activity in GRS 1915+105. 

In LHS and IS which we consider in our study, only a small part of the disk emission 
component is seen directly. The energy spectrum is dominated by a Comptonization com- 
ponent presented by a power law. To calculate the total normalization of the "seed" disk 
blackbody component we model the spectrum with a Generic Comptonization model [BMC 



XSPEC model, see details in Titarchuk, Mastichiadis & Kylafis (1997)] which consistently 



convolves a disk blackbody with a Green function of the Compton Corona to produce the 
Comptonization component. We argue that the disk emission normalization calculated using 
this approach produces a more accurate correlation with respect to the correlation with the 
direct disk component which was obtained using the additive model, multicolor disk plus 



power law [see e.g. McClintock & Remillard (2006) 



This Paper is a continuation of the study of index-QPO and index-seed photon normal- 
ization correlations in BH sources started in ST07 and ST09. Particularly here we present 
a study of the index-seed photon normalization (disk flux) correlation observed from GRS 
1915+105 when it evolves from LHS to HSS. The description of RXTE data-set used is 
given in §[2] We have analyzed a broader sample of state transitions from GRS 1915+105 
and we found a diverse phenomenology for index evolution through a transition. In £|3] we 
provide a detailed description of state transitions analyzed in this study. In §4 we discuss 
and interpret the results of our observational study. Specifically in §4 we consider the effect 
of the bulk motion Comptonization in the inner part of the accretion flow on the index 
evolution during a state transition. Also we show that the index saturation effect is a direct 
consequence of the existence of this inner bulk motion region and, therefore, can be consid- 
ered as an observational signature of the converging flow (black hole). Furthermore in §4 we 



discuss the TF04 model and the Monte Carlo simulations by Laurent & Titarchuk (2009) 



(in preparation) in which the observable index evolution with rh has been already predicted. 



Conclusions follow in §[5} 



Observations and Data Reduction 



In the present Paper, we have used publicly available data of the RXTE observatory 
obtained from January 1997 to April 2006. In total, our study includes 107 observations 
made at different BH spectral states (LHS, IS, HSS) of the system. Data sets were selected to 
represent a complete rise-middle-decay track of bright X-ray activity episodes behavior along 
bright radio flaring events (Si^ghz >250 mJy). Therefore we have chosen powerful (>250 
ASM counts/s) flaring episodes of GRS 1915+105 with a good coverage of simultaneous 
radio/X-ray observation. In the past some of these data of spectral transitions in GRS 



1915+105 were analyzed by 


Trudolyubov et al. 


(1999) 


Trudolyubov 


(2001 


), 


Muno et al. 


(1999 


), 


Reig et al. 


(2000 


), ST07 and 


Rodriguez et al. 


(2008 


) for the 1997 - 


1998 and 2005 



- 2006 transitions respectively. 

Standard tasks of the LHEASOFT/FTOOLS 5.3 software package were utilized for 
data processing using methods recommended by RXTE Guest observer Facility according to 
"The RXTE Cook Book" (http://heasarc.gsfc.nasa.gov/docs/xte/recipes/cook_book.html). 
For spectral analysis we used PCA Standard 2 mode data, collected in 3 - 20 keV energy 
range. The standard dead time correction procedure has been applied to the data. To con- 
struct broad-band spectra, data from HEXTE detectors have also been used. We subtracted 
background corrected in off-source observations. Only HEXTE data in 20 - 150 keV energy 
range were used for the spectral analysis in order to exclude the channels with largest un- 
certainties. The HEXTE data have been re-normalized based on the PCA. The data are 
available through the GSFC public archive (http://heasarc.gsfc.nasa.gov). In Tables 1-6 we 
list groups of observations covering the complete dynamical range LHS-(IS)-HSS-(IS)-LHS 
of the source evolution during flaring events. We present here period ranges MJD=50462 - 
51081 and MJD=53382 - 53852, as different types (samples) of bright X-ray activity, with 
transitions between hard and soft states. Two selected data sets have different patterns of 
radio/X-ray behavior and of light curve shapes. 

We also use public GRS 1915+105 data from the All-Sky Monitor (ASM) on-board 



RXTE (Swank 1999). The ASM light curves (2-12 keV energy range ) were retrieved from 



the public RXTE /ASM archive at HEASARC ( |http : / /xte. mit . edu / ASM _lc .htmip . 

The monitoring Ryle Radio Telescope (15 GHz) data in the 1997 - 2006 period were 
kindly provided by Dr. Guy Pooley. The technical details of the radio telescope are described 



by Pooley & Fender (1997) 
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2.1. Spectral analysis 

2.1.1. BMC and iron line components of the model spectrum 

The broad-band source spectra were modeled in XSPEC with an additive model consist- 
ing of two BMC: a BMC with high energy cut-off (BMC1 component) and BMC2 component: 
wahs* (bmc + bmc* highecut). We also use a multiplicative wabs model to take into account 
of an absorption by neutral material. The wabs model parameter is an equivalent hydrogen 
column Nh- Systematic error of 1% has been applied to the analyzed X-ray spectra. 

The BMC model describes the outgoing spectrum as a convolution of the input "seed" 
blackbody-like spectrum, whose normalization is Nb mc and color temperature is kT, with 
the Comptonization Green's function. Similarly to the ordinary bbody XSPEC model, the 
normalization Nb mc is a ratio of the source (disk) luminosity to the square of the distance 



.10 39 erg/sA d 

The resulting spectrum is characterized by the parameter log(A) related to the Comptonized 
fraction / as / = A/(l + A) and a spectral index a — T — 1. There are several advantages 
of using the BMC model with respect to other common approaches used in studies of X-ray 
spectra of accreting compact objects, i.e. sum of blackbody/multi-color-disk and power- 
law/thermal Comptonization. First, the BMC, by the nature of the model, is applicable to 
the general case where there is an energy gain through not only thermal Comptonization but 



also via dynamic (bulk) motion Comptonization (see Shaposhnikov & Titarchuk 2006, for 
details). Second, with respect to the phenomenological powerlaw model, the BMC spectral 
shape has an appropriate low energy curvature, which is essential for a correct representation 
of the lower energy part of spectrum. Our experience with powerlaw components shows 
that the model fit with this component is often inconsistent with the Nh column values 
and produces an unphysical component "conspiracy" with the highecut part. Specifically, 
when a multiplicative component highecut is combined with BMC, the cutoff energies E cut 
are in the expected range of 20~30 keV, while in a combination with powerlaw, E cut often 
goes below 10 keV, resulting in unreasonably low values for photon index. As a result, 
the implementation of the phenomenological powerlaw model makes much harder, or even 
impossible to correctly identify the spectral state of the source, which is an imminent task 
for our study. Third, and even a more important property of the BMC model, it calculates 
consistently the normalization of the original "seed" component, which is expected to be a 
correct mass accretion rate indicator. Note the Comptonized fraction is properly evaluated 
by the BMC model. 

We consider a scenario related to our spectral model (see Fig. [T]) where the Compton 
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cloud along with converging flow are located in the innermost part of the source and the 
Keplerian disk is extended from the Compton cloud (CC) to the optical companion (see e.g. 



TF04). An iron K^-line (laor) component (Laor 1991) was included in our model spec 



trum. To summarize the spectral model parameters are the equivalent hydrogen absorption 
column density Nh] spectral indices a±, a 2 (photon index T = a + 1); color temperatures 
of the blackbody-like photon spectra kT\, kT 2 ; log(4.i), log(4 2 ) related to the Comptonized 
fractions / 1; / 2 [/ = -4/(1 + 4)]; normalizations of the blackbody-like components N bmc i, 
Nbmc2 for the BMC1 and BMC2 components of the resulting spectrum, respectively. 

We find that color temperatures kT\ and kT 2 are about 1 keV for all available data 
and thus we fix values of kT\ and kT 2 at 1 keV. An equivalent hydrogen absorption column 



density was fixed at the level of Njj = 5 x 10 22 cm -2 (Trudolyubov 2001). When the 
parameter log(4) > 1 we fix log(4) = 2 (see Tables 1-6), because the Comptonized fraction 
/ = 4/(1 + 4) — »1 and variation of 4 does not improve the fit quality any more. 

During LHS the BMC2 component is often very low or barely detectable [see Table 1 
(MJD=50462-50561), Table 3 (MJD=51067-51081), Table 6 (MJD=53829-53852) and Fig.g 
and panel "IS"]. This observational fact is in agreement with scenario of BH spectral state 
transition (TF04). During LHS the spectrum is characterized by a strong hard power-law 
component. In other words the energy spectrum is dominated by a Comptonized compo- 
nent seen as a power-law hard emission in energy range from ~ 10 to ~ 70 keV while the 
disk emission remains weak (LHS, IS), because only a small fraction of the disk emission 
component (1 — /) is directly seen (Fig. [2j panel "IS"). 

Although during IS (Fig. [2j panel "2S" ) the contributions of two BMC components to 
the overall spectrum are of the same order sometimes we can barely identify one of these 
components [see these cases in Table 2, (MJD=50737, 50743), Table 3 (MJD=50908, 50909, 
51003), Table 5 (MJD=53718) when N 6mc2 <N 6mc i]. On the other hand the model with 
two BMC components are really needed in the most cases of the intermediate state (IS) and 
high-soft state (HSS). In Figure [3] we demonstrate that the fit qualities are unacceptable 
when the only one BMC component is included in the spectrum. Specifically, for IS-HSS 
observation 91701-01-11-00 on 18 May 2005 the values of x 2 ed =12.3 for 75 d.o.f. (bottom left 
panel) for wabs*bmc*highecut model. However x 2 is significantly improved when the second 
BMC component is included in the model. For wabs * (bmc + bmc* highecut) Xred = 3-28 for 
76 d.o.f. (see the related count spectrum along with the model and the related residual in 
the central bottom panel of Fig. |3|. Ultimately we achieve a remarkable agreement with the 
data using wabs * (bmc + bmc * highecut + laor) model for which x 2 ed =1.04 for 73 d.o.f. In 
Figure [3] (top and right bottom panels) we show the data along with the best-fit spectra and 
their components for our two component BMC model (see Table 4 for the best-fit parameter 
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values). 

In HSS and very soft state (VSS) (MJD=53515-53640, Table 4) the soft luminosity is 
high and the spectrum is dominated by a thermal component (T2 >3.7). Note the hard 
power-law component is barely seen in VSS (see Fig. § panels "4S"). We find a broad 
emission line between 6 and 7 keV in the IS and HSS spectra. Up to the date the iron K a 
emission line in IS and HSS of GRS 1915+105 was detected with ASCA, Chandra, XMM- 



Newton, BeppoSAX by Kotani et al. (2000), Martocchia et al. (2002), Miller & Homan 



(2005 ). However the determination of the iron line shape with RXTE is a problem because of 
the low energy resolution. As a first trial we added a Gaussian component to fit the spectrum 
varying the width and normalization of the line. Fits using the Gaussian always produce 
residuals around 7 keV while fits with XSPEC Laor model do not have such a problem. 
Thus through the Paper we incorporate the Laor- line profile to fit the line component. 
The line feature has a statistical significance of (3-10) a depending on the spectral states. 
This line feature is variable with time average intensity of 2xl0 -11 erg/s/cm 2 and exhibits 
equivalent width (EW) in the range of 50-600 eV across the data. We found that adding 
the laor-line component significantly improves the fit quality of IS and HSS spectra. Fitting 
an IS spectum (e.g. for 90105-05-03-05 observation) without the line components leads to 
Xred=l-24. When the line component (Laor) is included the fit quality becomes much better, 
^=1.01. The fit of HSS spectrum (901050-08-02-00) without the iron-line component is 
unacceptable, x 2 red =2.^\, and Xred = l-24 when the line component is included. The best-fit 
parameters of the source spectrum and values of Xred including d.o.f are presented in Tables 
1-6. 



2.1.2. Observational evidence of "blackbody-like" component peaked at ~ 20 keV in eight 

IS spectra 

The adopted spectral model shows a very good performance for 99 cases among 107 
spectra used in our analysis. Namely, the value of reduced x 2 -statistic Xred = xV-^W, 
where Ndof is a number of degree of freedom for a fit, is less or around 1.0 for most obser- 
vations. However for 8 observations of intermediate state the fit of the data with the model 
wabs* (bmc*highecut+bmc+ laor) is not so good, Xred reaches 1.5 and even higher. We found 
that in the residual of data vs model there is a characteristic bump around 20 keV (see left 
bottom panel of Fig. [4]) which can be fitted by blackbody-like shape of color temperature 
about 4.5 keV (see right bottom and top panels of Fig. [1] and Table 7 for values of the 
best-fit parameters). This "high-temperature BB" component is strong in each of the eight 
observations and its EW varies from 400 to 700 eV. We discuss a possible origin of this "BB" 
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component in §4. 



2.2. Timing analysis 

The RXTE light curves were analyzed using the powspec task. For the timing analysis 
in 2 - 30 keV energy range, we use the RXTE /PC A data in the binned and event modes 
containing X-ray events below 13/15 keV and above 13/15 keV for 1997/2005 data sets 
respectively. Specifically, depending on RXTE epoch, the channel ranges (0-35) for binned 
and (36-255) for event modes relate to energy bands 1.94-12.99 keV (binned) and 13.36- 
103.19 keV (event) for epoch 3 (1997-1998 data set), and relate to 2.06-14.76 keV (binned) 
and 15.18-117.86 keV (event) for epoch 5 (2005-2006 data set). 

The time resolutions for event and binned modes are 1.52 x 10~ 5 s and 8 x 10~ 3 s, 
respectively. The observational exposition time periods vary from 1.5 to 10 ks. Thus for all 
of these observations we can obtain power spectra in the wide frequency range (roughly from 
0.001 Hz up to 100/10000 Hz for binned and event modes respectively). These frequency 
ranges allow us to produce power spectra for all studied cases in the 0.1-100 Hz frequency 
range. We subtracted the contribution due to Poissonian statistics and Very Large Event 
Window for all power density spectra (PDS). 

The data analysis of the PDSs was performed using a simplified version of the diffusion 



model [see Titarchuk, Shaposhnikov & Arefiev (2007), hereafter TSA07, and Titarchuk & 



Shaposhnikov (2008)] in which the PDS continuum shape at frequencies below the driving 
frequency can be approximated by empirical model Px ~ (1.0+(x/x*) 2 )~~ in (KING model in 
QPD/PLT). Following TSA07, the break frequency found in the PDS is related to a diffusion 
time of perturbation propagation while the QPO low frequency is an eigenfrequency of the 
volume (magnetoacoustic) oscillation of the medium (in our case it is a Compton cloud). 
Note, TSA07 demonstrated that the diffusion model as a linear superposition of Lorentzians 
related to the eigenvalues of the diffusion problem can be also presented by the continuos 
shape which is flat below break frequency and power law at frequencies above the break. 
Given these asymptotic forms of PDS at low and high frequency limits they named their 
diffusion model as a white-red noise (WRN) model. For the quasi-uniform perturbation 
source distribution the slope of the PDS power-law part depends on the law of viscosity in 
the corona or in the disk (see details in TSA07). The parameters of this WRN diffusion 
model are the break (diffusion) frequency and the index of the power-law distribution of the 
viscosity over the radius. To fit the QPO features, we use Lorentzian shape. We quote the 
Lorentzian centroid as a QPO frequency. 
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3. Observational results 
3.1. Evolution of spectral properties during state transitions 



Observations of Galactic BH X-ray binaries reveal diverse spectral and dynamic phe- 
nomenologies. The evolution of a BH binary is usually described in terms of spectral states. 
There are several flavors of BH state definitions in literature, which slightly differ in BH 



state definitions and terminology (see, for example 


Hemillard & McClintock 


2006 


Belloni 


et al. 


2000 


2005 


Klein- Wolt & van der Klis 


200S 


). To distinguish different states, the 



properties observed in the energy spectrum and Fourier power density spectrum (PDS) are 
utilized. As we have already emphasized in the introduction section we use, in our study, 
the general state classification for four major BH states: low-hard (LHS), intermediate (IS), 
high-soft (HSS) and very soft state (VSS). 



The general picture of LHS-IS-HSS transition is illustrated in Figure [5] where we bring 
together spectra of LHS, IS, HSS and VSS to demonstrate the source spectral evolution 
from low-hard to soft states. We should emphasize different shapes of the spectra for the 
different spectral states. In the LHS spectrum the Comptonization component is dominant 
and the blackbody (BB) component is barely seen in 3-150 keV energy range. The IS and 
HSS spectra are characterized by a strong soft BB component and a power law extended up 
to 150 keV. In VSS the soft BB component is dominant and the power- law component is 
relatively weak with respect of this in IS and HSS. 

In the RXTE data of GRS 1915+105 observations there are long periods when the pho- 
ton index T\ and normalization N^mci of the hard BMC monotonically increase (or decrease) 
with time. We call these periods as long transition periods. The days, when the source 
X-tay flux starts to increase while it is still in the LHS, can be considered as a beginning 
of the rise transition. In these times the energy spectrum is characterized by low index 
values Ti ~1.7 and the thermal component is at low level or undetectable at all. In Figure 
[6] from top to bottom we show an evolution of flux density Si^ghz at 15 GHz (Ryle Tele- 
scope), RXTE/ ASM count rate, BMC normalization and photon index V during the 1997 
rise transition of GRS 1915+105 (MJD 50500 - 50700). Red/black points {for two low pan- 
els) correspond to hard/soft components with Y\ and T2, respectively. In the bottom we plot 
the photon index V versus the BMC normalization (left) and Comptonized fraction (right) 
for this transition. Here red triangles/black circles correspond to hard/soft components with 
Ti and T 2 , correspondingly. One can see that in the beginning of this transition the result- 
ing spectrum consists of one Comptonization component whose photon index Ti steadily 
increases from 1.7 towards the softer states and finally saturates at the value of 3. The 
Comptonization fraction / = A/(l + A) of the hard component, related to index Ti, shows 
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a sign of decreasing towards the softer state. When the RXTE /ASM count rate exceeds 50 
counts/s the soft Comptonized component appears with a weight which is comparable to 
that of the hard component. The photon index of the soft component T2 saturates to the 
level of 4.2 when the BMC normalization (disk flux) increases. The Comptonization fraction 
/ of the soft component is about 0.5 and higher. 

As seen from Fig. [6] the start of this rise transition coincides with active phase of X-rays, 
and of radio emissions which exceed 10 ASM counts/s and 50 mJy levels respectively in the 
1997 rise transition. Around MJD 50580 day the source reaches the HSS (when Ti ~ 3). 
Then a long HSS period from MJD 50580 to 50700, when Y\ stays almost the same, is 
followed by the state transition to IS during which Yi decreases to 2.5 (see Fig. [7]l8l. 



We see a similar behavior of X-ray, radio fluxes and photon indices during the 2005 
bright X-ray episode. The only difference of that with the 1997 episode was that the 2005 rise 
started at the intermediate state and went very quickly to HSS (see Fig. g[l0|). After MJD 



53800 the source came back to IS-LHS when r x ^2 (see Fig. 10 and Table 6). Note typical 



X-ray and radio fluxes during IS are about 40-60 ASM counts/s and <50 mJy respectively. 



3.2. Observational (correlated and non-correlated) characteristics of X-ray 

and radio emissions 



In fact, we do not find any correlation of X-ray and radio fluxes when the source changes 
its spectral states. Also we do not find a correlation of radio activity with the X-ray photon 



index (see Fig. 11). However we find a strong correlation of the iron line EW with radio 
flux density Sisghz at 15 GHz (see Fig. 12). In Figure 12 we also include points which 



have been recently reported by Neilsen & Lee (2009) who have analyzed archival HETGS 



(High Energy Transmission Grating Spectrometer) observations of GRS 1915+105 from the 
Chandra X-ray Observatory. 

The prominent HSS events were observed during 2005 - 2006 observations around MJD 
53490 and MJD 53690 (see Figs. @[l0|). The 2005 - 2006 observations confirm the index 
evolution vs BMC normalization (disk flux) found in the 1997 - 1998. The indices of the 
hard and soft components indeed increase and then saturate at values of ^3 and 4.2 respec- 



tively (see Fig. 10). Index Y\ started a saturation at lower values of BMC1 normalization 



(presumably proportional to disk mass accretion rates) than that were in 1997. 

It is also worth noting a so called "pivoting" effect, i.e. when inequality N bmcl > N bmc2 
switches to N bmc2 > N bmcl or vice versa, which is seen in the 1997-1998 and 2005 - 2006 
observations. One can see this composite pivoting picture combining Figs. [6 10 In fact, 
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these pivoting points correspond to the spectral transitions between adjacent states LHS-IS 
to HSS and vice versa. 



In Figure 13 we collect all data points for the index-normalization correlation for rise 
and decay stages. We do not find much differences in the correlation patterns related to the 
rise and decay transitions (compare left and right panels) in contrast that ST09 found in 
other BHs. 

We also find that the photon index of X-ray spectrum is tightly correlated with the 



quasi-periodic oscillations (QPO) frequency (see Fig. 14) which can be considered as a 
strong argument that QPOs and X-ray Comptonization spectrum emerge from the same 
geometrical configuration (Compton cloud). However the flux density S\ 5 ghz and QPO 
frequency are not correlated with each other when the source changes its spectral states. 



In Figure 15 we show an evolution of the flux density Si^ghz at 15 GHz (Ryle Telescope), 
RXTE I ASM count rate and vqpo during 1997 (left column) and 2005 (right column) rise 
transitions. The left column panel demonstrates the presence of QPO during a low radio 
activity (<30 mJy). The right column panel shows an example of the presence of QPO when 
the radio flux is high (~100 - 200 mJy). Given that the quasi-periodic oscillations (QPOs) 
of X-ray flux are present independently of the radio flux level, we can conclude that the 
radio appearances and QPO phenomenon are not closely related and probably the radio and 
X-ray (oscillating) emission areas have different origins in the source. 



3.3. Evolutions of energy and power spectra during a minor X-ray/radio 

flares 



In Figure 16 we show the details of a typical evolution of X-ray timing and spectral 



characteristics for minor X-ray/radio flares. In the top panels of Figure 16 we show the 
flux density S^ghz at 15 GHz (Ryle Telescope) and the RXTE / ASM count rate during the 
2005 rise transition stage (see also Fig. |9|. Red points A, B and C on the panel of the 
RXTE /ASM count rate vs time correspond to moments at MJD=53416, 53422 and 53442 
(before, during and after the minor X-ray/radio flare) respectively. Points A and C were 
chosen as the nearest possible points to point B (taking into account the time-table of archive 
data). Point B corresponds to the maximum of radio flux of 300 mJy and EW of 600 eV. 
Note that QPO centroid frequency before the flare (at point A) is at 1.8 Hz and shifts to 0.9 
Hz (point C) after the flare. 

PDSs (left bottom column) are plotted versus the energy spectrum (right bottom col- 
umn) for three points A (top), B (middle) and C (bottom) of the X-ray light curve. There 
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are QPOs at A and C points (Al, CI panels) but there is none at B point (Bl panel), at 
the X-ray flare peak. For the photon spectra (right bottom column) red points stand for 
observational data, while the model is shown by components with blue line for BMC1, black 
line for BMC2 and dashed purple line for the laor line component. Note that the spectral 
characteristics undergo noticeable changes during X-ray/radio flare. Specifically at the flare 
peak (point B) the total flux increases at least by a factor of 1.5 with respect to that before 
the flare. Although photon index of BMC1 component T\ changes from 2.9 (A and C points) 
to 3.0 (B point) respectively. 

We also studied the energy dependence of the PDS shape and integrated power vari- 



ability as a function of the photon energy. In Fig. 16 {left bottom column pannel) we show 
two power spectra for two energy bands 2-15 keV (red) and 15-30 keV (blue). One can see 
that PDSs weakly depend on the energy band. In particular, a value of the low frequency 
QPO vqpo is the same for the low energy and high energy PDSs. 



4. Interpretation and discussion of observational results 

Before to proceed with the interpretation of the observations let us to briefly summarize 
them as follows: i. The spectral data of GRS 1915+105 are well fit by two (soft and hard) 
BMC components for most of analyzed IS and HSS spectra (see Fig. |3| while LHS spectra 
essentially require only one BMC component, the soft BMC component is very weak (see 
Tables 1, 3-4, 6 and panel SI in Fig. [2]). ii. In addition to two BMC components 8 IS- 
HSS spectra require an extra component which can be fitted by " high temperature BB-like" 
profile (see Fig. [4] and Table 7). iii. The Green's function indices of each of these components 
rise and saturate with an increase of the BMC normalization (disk flux). The photon index 
saturation levels of the soft and hard components are about 4.2 and 3 respectively (see Fig. 



13). iv. We also find a tight positive correlation of QPO frequencies with the index (see 



Fig. 14) and consequently that with the disk flux. vi. The iron line EW correlates with 



the radio flux (see Fig. 12). vii. QPO appearances and their frequency values are not 



correlated with radio flux when the source undergoes the spectral changes from IS to HSS 



(see Fig. 15). vii. We also do not find any correlation between X-ray and radio fluxes and 
X-ray power-law index (see Fig. [TTj ). viii. Although we find changes of power and energy 
spectra during a minor X-ray/radio flares when QPO features disappear in PDS and energy 



spectrum becomes softer than that was before and after the flare (see Fig. 16) 



-14- 



4.1. Index-QPO and index-m correlations. Index saturation 

We confirm the index-QPO correlation in GRS 1915+105 previously found by V03 and 



ST07. This correlation was indeed predicted by Titarchuk, Lapidus & Muslimov (1998), 
hereafter TLM98, who argued that the transition layer [Compton cloud (CC)] formed be- 
tween the Keplerian disk and the central object (NS or BH), contracts and becomes cooler 
when the disk mass accretion rate rh increases. The observational effect of the CC contrac- 
tion were later demonstrated by ST06, TSA07, TS08 and Montanari, Titarchuk fc Frontera" 



| ( |2009 ) in Cyg X-l and XTE J1650-500 respectively. As a result of the transition layer 



(TL) contraction the QPO low frequency u^, which is presumably the TL's normal mode 
oscillation frequency, increases with m given that vl is inversely proportional to the TL size. 
On the other hand the index monotonically increases when the TL (CC) cools down. TF04 
provided the details of the index-QPO correlation model where they pointed out that this 
correlation is a natural consequence of the spectral state transition. 

In this Paper we have firmly established the index correlation with ul along with the 
index saturation vs the BMC normalization Nb mc (Eq. [T]) for the soft and hard Comptonized 
components of the X-ray spectra of GRS 1915+105 (see Fig. 13). Below we show that N bmc 



is actually proportional to mass accretion rate in the disk. Namely the disk flux L (as a 
source of soft photons for Comptonization, see e.g. Fig. [T] for the geometry of soft photon 
illumination of Comptonized region) can be represented as 

GM bh M 

L = - = r){r*)m d L Ed . (2) 

Here R* = r*Rs is an effective radius where the main energy release takes place in the disk, 
Rs = 2GM/c 2 is the Schwarzschild radius, r\ = 1/(2?%,), rhd = Md/M crit is dimensionless 
mass accretion rate in units of the critical mass accretion rate M crit = Lgd/c 2 and L^d is the 
Eddington luminosity. 



On the other hand 



ATiGMm p c 



<7t 

i.e. LEd oc M and thus using Eqs. ppl) we obtain that 



(3) 



L oc r](r*)mdrn. (4) 



In HSS when the inner disk radius R* reaches its lowest value R*>3Rs, the efficiency of 
the gravitational energy release 77(7%) reaches its highest value and thus the disk flux increases 
only when the disk mass accretion rate increases (see Eq. El). Given that BMC normalization 
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Nbmc is proportional to rhd in HSS the observational effect of the index saturation with iVb mc 
is translated to the index saturation with rhd- 

First we interpret the index saturation related to the hard Comptonization component 
(BMC1). We suggest that this BMC1 component of the emergent spectrum is presumably 
originated in the converging flow onto a compact object, in our case to the BH (see Fig. [I]). 
In fact, in HSS the plasma temperature of the accretion flow is comparable with the color 
temperature of the disk photons (see TF04). Thus, in order to explain the high energy tail 
observed in HSS of BH sources, one should assume either an unknown source of high energy 



non-thermal electrons [see e.g. Coppi (1999)] or consider effects of energy transfer from the 



converging flow electrons to the photons emitted from the innermost part of the accretion 
flow. 

Optical depth of the converging flow r is proportional to rhd if one assumes that disk 
accretion flow continuously goes to the converging flow and there are no other components 



in the accretion flow [see e.g. a model of two component accretion flow by Chakrabarti & 



Titarchuk (1995), hereafter CT95]. This effect of the index saturation vs optical depth 



of the bulk flow (BM) r was first predicted by Titarchuk & Zannias (1998) and then it 



was subsequently reproduced in Monte-Carlo simulations by Laurent & Titarchuk (1999), 
hereafter LT99. 

It is worth noting that the index saturation effect is an intrinsic property of the bulk 
motion onto a BH given that the spectral index a = T — 1 is a reciprocal of the Comptoniza- 



tion parameter Y [see this proof in ST09 and Bradshaw et al. (2007)] which saturates when 
the BM optical depth, or M, increases. In fact, the Y-parameter is a product of the average 
photon energy exchange per scattering 77 and the mean number of photon scattering N sc , i.e. 



Y = T)N SC . For the thermal Comptonization case, Y 
f] = AkT/m e c 2 and N sc ~ r 2 for r ^> 1 (see e.g. 
thermal Comptonization spectral index is 

a ~ [(4kT/m e c 2 )T 2 



(AkT / m e c 2 )r 2 given that in this case 



Rybicki & Lightman 1979) and, thus, the 



(5) 



In the case of converging flow, the preferable direction for upscattered photons is the direction 
of bulk motion onto the BH, i.e along the radius. Note that the fractional photon energy 
change is 

AE/E = (1 - MVc)/(l - fJL 2 V R /c). 

where ii\ and 112 are the cosines of the angles between the direction of the electron velocity 
n = V^/Vr and direction of incoming and outcoming (scattered) photons respectively. 

The number of scatterings of the up-Comptonized photons N sc can be estimated as a 
ratio of the radial characteristic size of the converging flow L and the free path I in the 
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direction of motion, namely N sc tx L/l — r given that AE/E has a maximum at {12 — 1 
for given fii and V#. On the other hand the efficiency per scattering for bulk motion flow 



r] oc 1/r when r ^> 1 Laurent & Titarchuk (2007), hereafter LT07] hence for bulk motion 



Comptonization, the Y-parameter does not depend on r for high values of r or dimensionless 
mass accretion rate rh. Thus one can conclude that the Comptonization parameter Y = r]N sc 
and hence the energy index a = Y~ l saturate to a constant value when optical depth (or mass 
accretion rate) of the BM flow increases. 

However the index saturation value is determined by the plasma temperature during 
a transition [see LT99]. The plasma temperature strongly depends on the mass accretion 
rate in the bulk motion region M& m and its illumination by the disk photons L (see TLM98 
and TF04). For higher M& m and L the plasma temperature is lower. The level of the index 
saturation decreases when the plasma temperature in the bulk motion increases (TF04). 
Thus the index saturation levels can be different from source to source depending on the 



strength of the disk. Looking at Figure 13 one can also notice that the index Ti starts 



its saturation at different values of BMC normalization (oc rhd) for different types of active 
episodes. In fact, the index should saturate with mass accretion rate in the converging flow 
fnbm which is a sum of the disk mass accretion rate rhd and that in sub-Keplerian flow 
(CT95). Hence one can argue that this lower value of rhd at which the index saturates can 
be a sign of the presence of extra (sub-Keplerian) component in the accretion flow onto BH 
in GRS 1915+105. 



Laurent & Titarchuk (2009), hereafter LT09, study the index-m correlation and also a 
modification of the disk blackbody spectrum due to Comptonization in the optically thick 
CC, which is formed due to accumulation of accretion matter in the TL. They indeed explain 
the saturations of the indices of the soft and hard components of the resulting spectrum. 
Specifically LT09 show that gravitational energy of the accretion flow is released in the 
optically thick and relatively cold TL when mass accretion rate rhd is higher than 1. The 
level of the index saturation depends on the radial velocity in the transition layer (TL). LT09 
also show that the observable saturation index of the soft BMC component T 2 ~ 4.2 can be 
reproduced in their Monte Carlo simulations for values of the TL radial velocity ^ 0.05 c. 



4.2. Physical origin of "high temperature BB-like" component"? 

In 8 of IS-HSS spectra we find an observational evidence of the bump around 20 keV 
which can be fitted by "~ 4.5 keV BB-like" profile (see Fig. [4] and Table 7). One can 
argue that this observable bump at 20 keV is a signature of the Compton reflection bump 



see e.g. Basko et al. (1974), ST80, CT95, and Magdziarz & Zdziarski (1995)]. But this 
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interpretation encounters difficulties given that the hard power-law tails of these spectra are 
too steep to form the Compton bump. Indeed, ST80 and later LT07 demonstrated that the 
Compton bump as a result of a photon accumulation due to downscattering of hard photons 
in the cold medium (for example disk) cannot be produced if the photon index of the incident 
hard photon spectrum T > 2. In fact, as one can see from Table 7, that in all spectra where 
we detect this ~ 20 keV feature the index of the hard BMC component r x > 2 (the indices 
vary between 2.5 and 3). 

In principle, this bump may also be a result of photoelectric absorption of the photons 
below 10 keV in the cold medium (disk) even if the incident spectrum is very steep. The 
photoelectric absorption cross-section a vh ~ (7.8 keV/E) 3 o"T, where E is photon energy 



and <tt is Thomson cross-section (e.g. CT95). However, Laming & Titarchuk (2004) and 



Rozahska & Madej (2008) show that the ionization of such a disk by the intensive X-ray 
radiation during IS-HSS invalidates the basic assumptions of the presence of the cold material 
in the innermost part of the source. Note the hard tail of X-ray spectrum in IS-HSS is formed 
in the converging flow (CF), i.e. in the innermost part of the accretion flow, because we see 
the CF signature (the index saturation) when the source goes to IS-HSS. 



On the other hand Titarchuk (2002) argued that the specific spectral and timing fea- 
tures of X-ray radiation could be seen in BH sources only. Particularly, he stated that the 
photon-photon interaction of the effectively upscattered photons results in the powerful pair 
production near a BH horizon. Indeed, a large fraction of the upscatterred photons going 
inward are deflected by the relativistic free-fall electron in the outward direction [the aber- 



ration effect of light, see e.g. Rybicki & Lightman (1979) and Appendix A]. These diverted 



upscattered photons of energy E up interact with incoming photons of energy E in flux and 
ultimately this interaction leads to the pair creation if the condition E up E in > (m e c 2 ) 2 is sat- 
isfied. Note that free-fall bulk motion with Lorentz factor 7 ^> 1 should be very close to 
horizon, i.e. 

R s 7 2 

where AR = R — R$ is radial distance to horizon. Thus the created positrons extensively 
interact with accreting electrons there and therefore the annihilation line photons are created 
and distributed over the relatively narrow shell near BH horizon. Specifically 

10\ 2 m 



AiTSxlO^-j - cm. (7) 

The proper energy (in the comoving flow frame) of annihilation line photons E$n should be 
seen by the Earth observer (in the zero frame) at the redshifted energy 



E = (l-R s /R)^E 5n - E ™ 



7 
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where (1 - R s /R) 1/2 « (Ai?/_R S ) 1/2 = l/ 7 (see according Eq. 

In other words the line energy displacement due to gravity as viewed by a far away 
observer in free space is z + 1 = 1/(1 — Rs/R) 1 ^ 2 = 7. 

A significant fraction of these annihilation line photons strongly gravitationally red- 
shifted can be directly seen by the Earth observer as a bump located at ~ 20 keV and 
presumably related to the representative value of 7 ~ 20. Laurent & Titarchuk (2009, in 
preparation) made extensive Monte Carlo simulations of the X-ray spectral formation in the 
converging flow taking into account photon-electron, photon-photon and pair-electron inter- 
actions. These simulations confirm our expectations that in some cases the emergent spectra 
of IS and HSS consists of the redshifted annihilation line located at ~ 20 keV, which can be 
fitted by "high temperature BB-like" profile, and also the simulated spectra are extended 



to energies of order of a few MeV [see Fig. [4] and Grove et al. (1998) for details of IS-HSS 
spectral components]. 



4.3. Radio— X-ray connection 



Migliari et al. (2005), hereafter MFK05 reported on correlations between radio lumi- 
nosity and X-ray timing (QPO) features in X-ray binary systems containing low magnetic 
field neutron stars and black holes. The MFK05 conclusions on the radio-QPO correlation 
based on observations of seven neutron star and one black hole GX 339-4. For GX 339-4 
they used data only in low-hard state before and after outburst. |Titarchuk fc Shaposhnikov 



[] (2005) based on the analysis of RXTE data from NS 4U1728-34 confirmed a correlation of 



X-ray and radio emissions with LF QPOs through all spectral states for this particular NS 
source. However, we do not find a real low-hard state in GRS 1915+105 for which photon 
index Ti should be about 1.5, as that in GX 339-4 (see ST09), and therefore we cannot 
confirm or refute the radio-QPO correlation in GRS 1915+105 LHS similar to that found by 
MFK05 in GX 339-4. 

On the other hand we find that LFQPOs do not correlate with the radio flux while 



they correlate with the hard component photon index Ti through IS and HSS (see Fig. 15 



and Fig. 14 respectively). Absence of correlation between radio luminosity and QPO can 
be explained by the different origins of these quantities. While the QPO phenomenon is 
probably related to the transition layer oscillations (see §4.1), it is confirmed by the index- 
QPO correlation, the radio emission is presumably originated in the wind or wide open 
jet. Furthermore, because the radio flux and iron line EW are strongly correlated (see 



Fig. 12) one can conclude that the iron line is also formed in the wind [see more on the 
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( 


2009 


) and 


Titarchuk, 



Shaposhnikov & Laurent (2009)]. It is also worth noting that the X-ray flux, and photon 



index do not correlate with the radio flux (see Fig. 11 ). It can be explained by the different 
mechanisms of the energy releases in X-ray and radio. X-ray radiation is presumably formed 
in the innermost part of the disk and in the transition layer (or Compton cloud) while the 
radio emission is presumably formed in the jet or winds which are probably launched at 



the outskirts of the disk [see e.g. Meyer-Homeister & Meyer (1999), Meyer et al. (2000), 



Meyer-Homeister & Meyer (2003) and TSA07] 



Using the aforementioned correlations of EW with radio flux we can suggest that some 
fraction of accretion flow may go to the outflow. Probably the powerful outflow is launched at 
outer parts of the accretion disk and it is not by chance that we can see this strong correlation 
of the iron line EW with the radio flux. Thus there are two ways for the matter to proceed: 
i. in the outflow if the local mass accretion rate in the disk exceeds the critical value (which 
is proportional to radius, see TSA07), ii. in the disk where the matter proceeds and it 
ultimately converges onto BH. This final stage of the accretion we observe as a saturation 
of index with mass accretion flow (converging flow signature). 

Whereas we do not find any correlation between radio and X-ray fluxes during the 
spectral transition, although, probably, we see an indication of X-ray-radio connection during 



a minor flare event. In Figure 16 we show the spectral and timing properties of X-ray emission 
of a typical minor X-ray/radio flare As one can see from this Figure (see low panels there) 
the PDS and energy spectra are different at the peak of the flare from those before and after 
the flare. Specifically the soft component is more pronounced but QPO features are not 
seen at the peak of the flare while they are present before and after the flare. Moreover the 
flat part (white noise) of the peak PDS is extended to higher frequencies (break frequency 
Uf, ~ 10 Hz at the peak vs i/& ~ 2 Hz before and after the flare). 



In §3.3| we reported the results of the study the energy dependence of the PDS shape 
and integrated power variability as a function of the photon energy. In Figure [16] we show 
two power spectra for two energy bands 2-15 keV (red) and 15-30 keV (blue). One can see 
that PDSs weakly depend on the energy band. We can suggest that the sizes of the photon 



emission areas Lcc related to these two energy bands are also the same. In fact, Titarchuk 



et al. 



(2001) argue that vqpo is proportional to the ratio of magneto-acoustic (plasma) 
velocity Vma and Compton cloud size Lcc an d hence one can conclude that the emission 
areas are the same because vq PO are the same for these two energy bands. 



Lightman & Eardley (1974), hereafter LE74, suggest that the thin disk is always un- 



stable in the inner region when radiation pressure dominates gas pressure. Using numerical 



simulations Lightman (1974) found that the innermost region of a disk around a BH is 
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secular unstable against clumping of the gas into rings which observational appearances can 
be seen by the Earth observer as X-ray-radio flares. We can speculate that an increase of 
the soft and hard components in the observable X-ray spectrum at the flare peak (compare 



panels B2 and A2, C2 of Fig. 16) can be a sign that the radiation pressure gets to dominate 
in the inner disk region. On the other hand the sign of the destruction of some part of the 
innermost part of the disk, as an effect of the high pressure instability, should be seen in 
the PDS. TSA07 argue that the break frequency u h in PDS is proportional to the diffusion 
frequency v& = l/t V i SC ~ v/R 2 where v is a viscosity, t V i SC is a viscous timescale and R is a 
radial size of the innermost part of the disk. Given that increases at the peak with respect 



of that before and after the flare (compare panels Bl and Al, CI of Fig. 16) it can imply 
that the size R decreases when increases, i.e. some part of the innermost part of the disk 
is probably destroyed. 

In terms of the diffusion theory, this disk instability arises in the inner region where 
the viscous stress W is a decreasing function of the surface density E and thus an effective 
diffusion coefficient of the nonlinear equation for £ becomes negative there [see Eqs. 4-5 
in LE74]. On the other hand Makeev & Titarchuk (2009, in preparation), hereafter MT09, 
obtain this disk instability as a solution the linear diffusion equation for £ and they do not 
specify any (ad hoc) assumption about the nature of the disk viscosity (cf. LE74). They just 
assume the power-law viscosity distribution over the disk and they use the TLM98 angular 
velocity distribution in the transition layer (TL). 

MT09 study the TSA07 model of the PDS formation and they find that this model 
predicts the existence of the two distinct zones within the TL of the black hole or neutron star, 
with oppositely different types of diffusion of perturbation taking place in each zone. A simple 
fact of the change of sign of the angular velocity derivative in the linear diffusion equation 
for E at the critical radius R ma x results in a negative diffusion coefficient of the equation in 
the interval R{ n < r < R ma x- Moreover the change of sign of the diffusion coefficient at R ma x 
leads to a turnover of the angular momentum transfer, changing it towards the central body 
of the accreting system, instead of being pushed outwards as in the Keplerian disk. One of 
the considered scenarios implies an unstable diffusion of the perturbations in the TL inner 
zone which might be an indication of the development of a X-ray flare followed by a radio 
flare. The absence of any quasiperiodic oscillations in the TL and raising Ub are possible 
indications of this instability. 
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5. Conclusions 

We concentrate our efforts on the study of the correlation between the spectral index 
and the accretion disk luminosity. We argue that the shape of the correlation pattern can 
contain the direct BH signature. Namely, we show both observationally and theoretically 
that the index saturates with mass accretion rate which is a signature of a converging flow. 
This index saturation effect can exist only in BH sources. Also this correlation pattern carries 
the most direct information on the BH mass and the distance to the source (see ST09). 

We compiled the state transition data from GRS 1915+105 collected with the RXTE 
mission. We examined the correlation between the photon index of the Comptonized spectral 
component, its normalization and the QPO frequency (see Figs.[§^|l3||l4|. 



The spectral data of GRS 1915+105 are well fitted by two (soft and hard) BMC com- 
ponents for most of analyzed IS and HSS spectra (see Fig. [3]) while LHS spectra essentially 
require only one BMC component. In addition to two BMC components 8 IS-HSS spectra 
require an extra component which can be fitted by " high temperature BB-like" profile. 
We suggest this "BB" component is probably a signature of the redshifted annihilation line 
formed in the very narrow shell near a BH horizon due to high photon compactness taking 
place during intermediate and high-soft states (see Fig. Eland Table 7). 



A remarkable result of our study is that the index - normalization (mass accretion 
rate) correlation seen in GRS 1915+105 is predicted by the theory of the converging flow. 
We demonstrate that a strong index saturation vs disk flux seen in the index-disk flux 



correlation (see Fig. 13) is an observational signature of the presence of the converging flow, 
which should only exist the BH sources. In other words, this index saturation effect provides 
a robust observational evidence for the presence of black hole in GRS 1915+105. 

We also find a tight positive correlation of QPO frequencies with the index (see Fig. 



14) and consequently that with the disk flux. 



Our comprehensive analysis of X-ray and radio emissions in GRS 1915+105 shows that 
QPOs are seen independently of radio activity of the source during the spectral transition 
from low-hard to high-soft state. Specifically these QPO features have been detected at any 
level of the radio flux and even when the radio emission is at the noise level (see left and 



right panels in Fig. 15 correspondingly). We also do not find any correlation between X-ray 



and radio fluxes and X-ray power- law index, (see Fig. 11). However, we establish a strong 
correlation between equivalent width of iron line and radio flux in binary GRS 1915+105 
(see Fig. ph. 
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A. The light aberration effect 



For a given electron velocity v moving inward along radius v = — v e r the directions of 
the ray in the observed laboratory and comoving electron frames n and n are related to 



aberration formula [see e.g. Rybicki & Lightman (1979) 

v fc — cos 9q 



cos 9 



(v j c) cos 6*o — 1 



(Al) 



where 7 = [1 — (v/c) 2 ] x l 2 is Lorentz factor, cos#o = — (n • e r ) and cos# = — (n • e r ) are for 
the laboratory (zero) and electron rest frames respectively. 

For 7 ^> 1, we can write that 1 — v/c — l/(2~f 2 ) and thus we can rewrite Eq. (Al) as 
follows 

-2sin 2 (V2) + l/(2 7 2 ) 



cos 9 



If # ~ &/7 then for b > 1 



2sin 2 (# /2) + cos# /(27 2 ) 



cos 9 ~ — < 



(A2) 



(A3) 



1 + b 2 

i.e. ray is directed against the electron direction. In fact, cos 9 ~ —1 for all 9q ^> I/7. On 
the other hand 

cos# = -(n- e r ) > (A4) 

if b < 1. 
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Converging inflow Bulk comptonization 



Fig. 1. — Schematic view of the proposed geometry for thermal and bulk Comptonization 
regions in the source hosting a BH with PL-like emission at high energies. The thermal plus 
bulk Comptonization spectrum (thermal plus bulk BMC1) arises in the innermost part of the 
transition layer (TL), where the disk BB-like seed photons are (thermally and dynamically) 
Comptonized by the in-falling material. Whereas the thermal Comptonization spectrum 
(thermal BMC2) originates in the outer part of the TL region. 
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Fig. 2. — Evolution of spectum shape of GRS 1915+105 during LHS, IS, HSS and VSS 
spectral states. Data are taken from RXTE observations 20402-01-11-00 (SI, V\ = 1.8, 
LHS), 91701-01-33-00 (S2, T 1 = 2.9, T 2 = 3.7, IS), 91701-01-11-00 (S3, T 1 = 2.7, T 2 = 4.1, 
HSS) and 91701-01-19-00 (SI T 2 = 4.2, VSS). Here data are denoted by red points, the 
spectral model presented with components are shown by blue, black and dashed purple lines 
for BMC1, BMC2 and laor components respectively. 
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Fig. 3. — The best-fit spectrum of GRS 1915+105 during HSS events of 2005 rise transition 
(observation 91701-01-11-00 ) in E * F(E) units (top) along with the count spectrum and 
Ax (bottom panels). Bottom panels: Left: the best-fit spectrum and Ax for the model fit, 
wabs*bmc*highecut, which include a single BMC component without an iron line component 
(xl ed =12.3 for 75 d.o.f.), center: the model fit, wabs* (bmc + bmc*highecut), which includes 
two BMC components without an iron line component (Xr ec 2 = 3-28 for 76 d.o.f.) and right: 
the same as the latter one but adding of the Laor line component, wabs * (bmc + bmc * 
highecut + laor), (x^ ed =1.04 for 73 d.o.f., see Table 4). Blue, black, purple lines stand for 
BMC1, BMC2 and Laor components respectively. 
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channel energy (keV) channel energy (keV) 

Fig. 4. — The best-fit spectrum during IS-LHS of 2005-2006 decay transition in EF(E) units 
(top) and in normalized counts units (bottom panel) with Ax for the 91701-49-00 observation. 
Left: fitting without modelling the high-temperature "bbody" component (xl ed =^-A?> for 73 
d.o.f.) and right: the best-fit spectrum and Ax, when the bump in residuals at ~ 20 keV is 
modelled by a "high-temperature bbody-like" component with x^ ed =1.06 for 70 d.o.f.. On top 
panel the data are denoted by red points, the spectral model presented with components are 
shown by blue, black, purple and green lines for BMC1, BMC2, laor and "high-temperature" 
bbody components respectively. 
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Fig. 5. — Four representative EF E spectral diagrams during LHS, IS, HSS and VSS spectral 
states of GRS 1915+105. Data are taken from RXTE observations 20402-01-11-00 (blue, 
LHS), 91701-01-33-00 (red, IS), 91701-01-11-00 (green, HSS) and 91701-01-19-00 (purple, 
VSS). 
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Fig. 6. — From Top to Bottom: Evolutions of the flux density S\^ghz at 15 GHz (Ryle Tele- 
scope), RXTE /ASM count rate, BMC normalization and photon index V in the beginning 
of the 1997 rise transition of GRS 1915+105 (MJD 50460-50700). Red/black points (for 
two last panels) correspond to hard/soft components with I\ and T 2 , respectively. Bottom: 
Spectral index V plotted versus BMC normalization (left) and Comptonized fraction (right) 
for this transition. Here the red triangles/black circles correspond to hard/soft components 
with Ti and 1^, correspondingly. Note that in most cases the normalization of the hard BMC 
component (BMC1) is higher than that of the soft component (BMC2) (see red points vs 
black points in BMC normalization-time panel). 
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Fig. 7. — From Top to Bottom: The same types of evolutions which are presented in Fig. 
[6] but now that for the middle of the 1997 rise transition of GRS 1915+105 (MJD 50650- 
50760). In most cases the normalization of the soft BMC component (BMC2) is higher than 
that of the hard component (BMC1) [see black points vs red points in BMC normalization - 
time panel and compare with that in Fig. |6). The pivoting effect occurs (i.e. N bmcl > N bmc2 
switches to N bmc2 > N bmc i around MJD 50690. 
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Fig. 8. — From Top to Bottom: The same types of evolutions which are presented in Figs. [6} 
[7] but now that for the 1997 - 1998 decay transition of GRS 1915+105. In most cases 
the normalization of the hard BMC component (BMC1) is higher than that of the soft 
component (BMC2) (see red points vs black points in BMC normalization-time panel). 
The pivoting effect occurs around MJD 50910. 
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Fig. 9. — From Top to Bottom: The same types of evolutions which are presented in Figs. 
[6]j8]but now that for the 2005 rise transition of GRS 1915+105. The pivoting effect occurs 
around MJD 53475. 
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Fig. 10.— From Top to Bottom: The same types of evolutions which are presented in 
Figs. [6]j9] but now that for the 2005 - 2006 middle and decay transitions. The pivoting effect 
occurs between MJD 53750 and MJD 53775. 
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Fig. 11. — BMC photon indices T (top), normalization (middle) and X-ray flux (bottom) 
vs flux density Si 5 ghz- Top and middle: Red and black points stand for hard and soft 
spectral components respectively for the 1997 - 1998 active episode. Crimson and blue 
points correspond to hard and soft spectral components respectively for the 2005 - 2006 
episode. Bottom: Blue circles and red triangles stand for 1997 - 1998 and 2005 - 2006. 
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Fig. 12. — The equivalent width (EW) of iron line in eV as a function of flux density S^qhz 
at 15 GHz (Ryle Telescope) in mJy for 1997/2005 data sets (blue/red points). Here we 



also include black points which have been recently found by Neilsen & Lee (2009) analyzing 



archival HETGS (High Energy Transmission Grating Spectrometer) observations of GRS 
1915+105 from the Chandra X-ray Observatory. 
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Fig. 13. — Photon index plotted versus BMC normalization for rise (left) and decay (right) 
transitions. Red/black points stand for hard/soft spectral components for the 1997 activity 
episode. Crimson and blue points correspond to hard and soft spectral components respec- 
tively for the 2005 episode. Points marked with rose oreol correspond to IS-HSS spectra 
fitted by the model which includes "high-temperature bbody" component (for details see 
Fig. |]and Table 7). 
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Fig. 14. — Photon index Ti (hard component) plotted versus QPO centroid for the 1997 rise 
transition from GRS 1915+105. 
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Fig. 15.- Evolution of flux density S 15GHz at 15 GHz (Ryle Telescope), RXTE / ASM 
count rate, vqpo during 1997 (left column) and 2005 (right column) rise transitions of 
GRS 1915+105. Here vq PO stands for the centroid frequency of the fundamental QPO. 
The left column panel demonstrates the presence of QPO even when for radio flux is low 
(<30 mJy). The right column panel shows an example of the presence of QPOs during 
prominent radio flux events (~100 - 200 mJy). These panels demonstrate that the QPO 
appearances are independent of the radio flux during LHS-HSS transitions. 
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Fig. 16. — Top: evolution of the flux density Si^ghz at 15 GHz (Ryle Telescope), 
RXTE /ASM count rate during the 2005 rise transition. Top lower panel: red points A, 
B and C mark moments at MJD=53416, 53422 and 53442 (before, during and after radio 
flares) respectively. Bottom: PDSs (left column) are plotted along with energy spectral dia- 
gram EF(E) (right column) for three points A (top), B (middle) and C (bottom) of the X-ray 
light curve. Power density spectra for the soft (red, 3 -15 keV) and hard (blue, 15 - 30 keV) 
energy bands are presented. There are QPOs at A and C points (Al, CI panels) but there 
is none at B point (Bl panel), at the X-ray flare peak. Here data are denoted by red points, 
the spectral model presented with components are shown by blue, black and dashed purple 
lines for BMC1, BMC2 and laor components respectively. 



